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The nonequilibrium dynamics of a binary Lennard- Jones mixture in a simple shear flow is investi- 
gated by means of molecular dynamics simulations. The range of temperature T investigated covers 
both the hquid, supercooled and glassy states, while the shear rate 7 covers both the linear and 
nonlinear regimes of rheology. The results can be interpreted in the context of a nonequilibrium, 
schematic mode-coupling theory developed recently, which makes the theory applicable to a wide 
range of soft glassy materials. The behavior of the viscosity ri{T, 7) is flrst investigated. In the 
nonlinear regime, strong shear-thinning is obtained, r/ ~ 7~"''^', with a{T) ~ 2/3 in the super- 
cooled regime. Scaling properties of the intermediate scattering functions are studied. Standard 
'mode-coupling properties' of factorization and time-superposition hold in this nonequilibrium situ- 
ation. The fluctuation-dissipation relation is violated in the shear flow in a way very similar to that 
predicted theoretically, allowing for the definition of an effective temperature Tcff for the slow modes 
of the fluid. Temperature and shear rate dependencies of T^h axe studied using density fluctuations 
as an observable. The observable dependence of T^a is also investigated. Many different observables 
are found to lead to the same value of TcS, suggesting several experimental procedures to access Toff. 
It is proposed that a tracer particle of large mass mtr may play the role of an 'effective thermometer'. 
When the Einstein frequency of the tracers becomes smaller than the inverse relaxation time of the 
fluid, a nonequilibrium equipartition theorem holds with (mtrvl^ — kBTcS, where Vz is the velocity 
in the direction transverse to the flow. This last result gives strong support to the thermodynamic 
interpretation of T^a and makes it experimentally accessible in a very direct way. 

PACS numbers: 64.70.Pf, 05.70.Ln, 83.60. Df 

L'immense thermometre refletait modestement la temperature du milieu ambiant. 
Vassili Axionov, L'ile de Crimee. 



I. INTRODUCTION 

Glasses are usually defined by the fact that their inter- 
nal relaxation time is larger than the experimental time 
scale. In simple molecular systems, the associated glass 
transition temperature corresponds to a very high vis- 
cosity, making it difficult to investigate experimentally 
the Theological properties of glassy systems. In complex 
fluids (e.g. colloids, emulsions) it is however possible to 
reach a glassy situation, in the sense of large relaxation 
times, with systems having viscosities or shear moduli 
that allow for a rheological investigationtl. Such materi- 
als have been described as 'soft glassy materials' in the 
literatureB. 

When it is quenched into its glassy state, a material 
is by definition out of equilibrium. An important fea- 
ture of this nonequilibrium situation is the absence of 
time translation invariance. Physical properties are a 
function of the time spent in the glassy phase, or wait- 
ing time tyj. This is best seen through the measurement 
of time-dependent correlations which depend both on i^, 
and on the time difference. These dependencies on 
are usually described as aging phenomenaEI. They have 
been jStudied experimentally and theoretically in great 
detain. Interestingly for the present work, the aging be- 
havior of several comnlepLJluids has recently been exper- 
imentally investigatecuLJ, showing interesting similari- 



ties with other, more standard, glassy systems. 

By submitting the system to a homogeneous, steady 
shear flow, a different kind of nonequilibrium situation 
is obtained. The steady shear flow, characterized by the 
shear rate 7, creates a nonequilibrium steady M,ate, in 
which time translation invariance is recoveredE3. The 
shear flow can therefore be used as a probe of the glassy 
system, with the convenient feature of having the shear 
rate 7 rather than the waiting time as a control pa- 
rameter. In fact, the inverse shear rate introduces a time 
scale in the problem, which plays a role similar to tw 
Interesting phenomena are thus expected to show up as 
soon as this new time scale competes with the relaxation 
time of the fluid. As a result, it is also interesting to 
study the 'supercooled' regime, which would correspond 
to an equilibrated situation in the absence of a shear 
flow. Moreover, this way of probing the nonequilibrium 
properties of glassy systems is possibly more relevant ex- 
perimentally than the aging approach, at least in the case 
of soft glassy materials. 

Recently, a general scenario was proposed for glassy 
systems subject to an external forcingE£l, based on the 
study of simple 'mean-field' models. The rationale of 
this approach is that the equilibrium dynamics of these 
models is equivalent to the 'schematic' mode-coupJiBC 
approach of slowing down in supercooled liquidsEjtJ. 
The study of their nonequilibrium dynamics can thus be 
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seen as-a nonequilihrium schematic mode-coupling ap- 
proactO. In this respect, it is interesting to note that 
the mode-coupJing theory was recently extended to the 
aging regimdla with results which . wfirc first predicted 
from the study of mean-field modclsll3c3. To our knowl- 
edge, the mode-coupling theory of supercooled fluid has 
not yet been extended to fluids under shear beyond the 
linear response of the supercooled regime, i.e. at equi- 
librium. However, in analogy to what was done for ag- 
ing or supercooled systems, it seems sensible to bypass 
this aspect and to carry out a direct comparison between 
the predictions from mean-field theories and experimen- 
tal or simulation results, with the hope that the general 
scenario is robust enough that the details of the system 
under study are relatively unimportant. 

Several earlier studies have been devoted, to simulating 
fluid under shear. Onuki and YamamotoEll concentrated 
their works on the understanding of the dynamics at the 
molecular level. Our study is focused on more global 
aspects, with the aim of providing same experimentally 
testable predictions. Liu and NageE3 also proposed to 
use the shear rate as a relevant coptrol parameter for 
jamming systems. Liu and coworkersES investigated some 
aspects of the nonequilihrium dynamics of a model sim- 
ilar to ours, with the difference of being athermal (zero 
temperature) in the absence of shear. Some of their re- 
sults are closely related to ours, as will be discussed in 
Sections |^ and VI . Mom-»henomenological approaches 
may be found in Refs.BBEJ. 

The aim of the present work is to check, on a realistic 
model of the fluid, some of the predictions that emerged 
from the earlier study of driven glassy systemalj. For 
completeness, the main results of this study will be briefly 
recalled in section In section III, we describe our mi- 
croscopic model for the fluid under shear. Sections IV 
and ^ describe our results for the dynamic properties 
of the fluid, both at the macroscopic (rheological) and 
at the microscopic (wavevector dependent correlations) 
scales. Sections ^ and |VIl| are devoted to one of the 
most important predictions of the mean-fleld scenario, 
namely the manner in which the equilibrium fluctuation- 
dissipation theorem is violated in such nonequilihrium 
systems and the resulting notion of an effective temper- 
ature. Section VIII proposes a numerical realization of 



a very simple experimental protocol to measure the ef- 
fective temperature through a nonequilihrium general- 
ization of the equipartition theorem. We discuss exper- 
imental and theoretical consequences of our results in 
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II. SUMMARY OF THE NONEQUILIBRIUM 
MODE-COUPLING RESULTS 

To make the paper self-contained, we briefly recall in 
this Section the main results and predictions obtained 
within the mean-field approach of Ref.E3. In this pa- 
per, a simple system, the p-spin mean-field model, was 



studied under conditions that involved a constant energy 
input from an external driving force, analogous to the 
effect of a shear flow on a fluid. The resulting station- 
ary nonequilihrium state was studied as a function of the 
temperature T and external drive, the intensity of which 
will be, for convenience, denoted by a. 

For (7 = 0, the model has an ideal glass transition at a 
finite temperature Tc where the relaxation time diverges. 
For T > Tc, the correlation functions are described by 
usual mode-coupling integro-differential equationsHj, and 
display the characteristic two-step relaxation predicted 
by these equations, the a- and /3-relaxations. Below Tc, 
the-system cannot equilibrate, and displays aging behav- 

iorka. 

Under a finite external drive, a ^ 0, the system is 
stationary at all temperatures. For T ^ Tc, and in the 
asymptotic limit cr ^ 0, the correlation functions re- 
tain the characteristic two-step shape of the equilibrium 
system, with an a-relaxation time that depends on a. 
For T < Tc, this relaxation time diverges as it —^ 0, 
while for T > Tc it goes to its equilibrium value in this 
limit. Although the driving force strongly influences the 
a-relaxation, it does so by keeping the shape of the decay 
of the correlation unchanged. This leads to the property 
that correlation functions may be collapsed by a simple 
rescaling of the time. This scaling property is analogous 
to the time-tcmpcrature superposition property found at 
equilibrium, and we shall denote this prediction by 'time- 
shear superposition property'. 

Based on a power dissipation argument, we were able 
for this simple system to define quantities equivalent 
to the viscosity and shear rate in a fluid under shear. 
This viscosity was found to have a characteristic shear- 
thinning behavior, decreasing for increasing shear rates 
as a power law. The shear-thinning exponent is equal to 
2/3 for T > Tc, while for T < Tc, the viscosity diverges 
as cr ^ 0, and the shear-thinning exponent appears to 
depend on temperature with a value between 2/3 and 1. 

The behavior of the effective temperature was also in- 
vestigated as a function of T and a. We recall here that 
the effective temperature Tcs is defined, in a systeia in- 
variant under time translation, by the relationshiplij 
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dC(t) 



fcBTeff(C) dt 



(1) 



Here R{t) and C{t) are, respectively, a response func- 
tion and the associated correlation function. In equilib- 
rium, the fluctuation-dissipation theorem can be written 
as Tes{C) — T. In the sheared system, the effective tem- 
perature was found, for T < Tc, to be a discontinuous 
function of C in the limit cr ^ 0. For C > q, one has 
Teff = T while for C < q, one flnds Teff(C) = T^g > T, 
denoting an effective temperature, larger than the bath 
temperature for long time scalealj. The parameter q is 
the plateau value of the correlation function: it is called 
the Edwards- Anderson parameter in the literature of dis- 
ordered systems, or nonergodicity parameter in the lan- 
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guage of mode-coupling theory. At finite shear, the ef- 
fective temperature also has an almost discontinuous be- 
havior as a function of C. For T > Tc, the discontinuity 
decreases and vanishes when a 0, so that equilibrium 
behaviour is recovered in this limit. 

An important outcome of theoretical studies of the 
concept of effective temperature for the slow modes of 
a glassy system, is that the value of the effective tem- 
perature is independent, at the mean-field level, of which 
observabl e |i s ch osen to compute correlation and response 
functionstSESnJ'EJ. This prediction, was recently chal- 
lenged within a simple trap modclEZI, with negative re- 
sults. We will see in Section VI that this important pre- 



diction is nicely verified in our realistic model, j— . 

The simple mean-field model studied in Ref.EJ does 
not allow any discussion of spatial dependency. However, 
characteristic features of the spatial dependence usually 
expected from the solution of mod e-cou pling like equa- 



tions will be investigated in Section VC 



III. 



MODEL AND DETAILS OF THE 
SIMULATION 



The system simulated in this work is a 80:20 mixture 
oi N = 2916 Lennard- Jones particles of types A and B, 
with interaction 



7al3 
''a/3 



(2) 



where a and /3 refer to the two species A and B. Interac.- 
tion parameters are chosen to prevent crystallization^. 
In all the paper, the length, energy and time units are the 
standard Lennard-Jones units uaa (particle diametec^, 
£AA (interaction energy), and tq = {mAcr\A/ saaY^'^^, 
where niA is the particle mass and the subscript A refers 
to the majority species. Particles have equal masses, and 
the interaction parameters are eab = 1-5 £aa, sbb = 
Q.^eaa, o-bb = 0.88 tJAA, <^AB = 0.8 aAA- With these in- 
teraction parameters between species, equilibriurrL-prop- 
erties of the system have been fully characterized^. At 
the reduced density p — 1.2, where all our simulations arc 
carried out, a 'computer glass transition' is found in the 
vicinity of Tc ~ 0.435. The slowing down of the dynam- 
ics, T > 0.45, is correctly described by mode-coupling 
theory, which breaks down when lowering further the 
temperature, T < 0.45Ea. The aging behavior of the 
system below this temperature has also been character- 
ized extensively,|-including the violation of the FDT in 
the glassy phaseE3. 

In order to study the system under a steady shear flow, 
the classical Newton equations are replaced by the so- 



called SLLOD equation; 
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- Pi ■ Vv, 



(3) 



where {pi,ri) are the momentum and position of parti- 
cle I, respectively. These eauations are integrated by a 
standard leapfrog algorithn£3, where time is discrctizcd. 
The value At = 0.01 was used throughout the simula- 
tions. Lees-Edwards boundary conditionsE3 are used in 
a cubic simulation box of linear size L = 13.4442. With 
these boundary conditions the flow is homogeneous and 
no instability, sucb-as^e.g. shear-banding, was observed 
in our simulationsOo. Note also that this is the shear 
rate 7 which is controlled in our simulations. The veloc- 
ity gradient is in the y direction, and the fluid velocity in 
the X direction, i.e. v = "fy e^- Constant temperature 
conditions are ensured by a simple velocity rescaling, of 
the z component of the velocities, at each time stepE3. 

The shear rate 7 naturally introduces a new time scale 
7~^ into the problem. Obviously, a simulation involv- 
ing a steady shear state is possible only if the available 
simulation time is significantly larger than 7^^. This lim- 
its our study to shear rates larger than typically IO'^tq, 
corresponding to 10^ time steps. The temperature range 
studied here is T e [0.15, 0.6], which covers the liquid and 
glassy phases, while the shear rates are 7 G [10"'', 0.1] 
which covers the linear and nonlinear response regimes 
of rheology. All the data presented in the following, even 
those below Tc, are obtained in a stationary state. This 
is unusual for systems below the glass transition tem- 
perature, which normally exhibit a non-stationary aging 
behavior. 



IV. MACROSCOPIC BEHAVIOR: NONLINEAR 
RHEOLOGY 

In this Section, we study the rheological behavior of 
the system. Our results are very similar to that obtained 
in many different soft condensed matter systems, mak- 
ing the Lennard-Jones mixture — which was originally 
designed as a model for simple metallic glasses — a rea- 
sonable one for more complex systems. 

Before discussing our results, let us briefly recall some 
possible behaviors for the flow curves (stress cr as a func- 
tion of shear rate) obtained in complex fluids. These 
resu^^fj^ are often represented in the phenomenological 
forr 
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(4) 



The case n = 1 corresponds to Bingham fluids, with a 
yield stress ctq, which vanishes for Newtonian systems. 
For|-p < 1, Eq. (|^) is often called the Hershel-Bulkeley 
lawtl. For zero yield stress (ctq = 0), Eq. (0) describes a 
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'power-law fluid', where the viscosity r] decreases with 
increasing the shear rate as a power law, rj ~ 7~", 
a = 1 — n. Although these type of flow curves are found 
in many different complex systems, they are not usually 
supported by a theoretical basis, and should therefore be 
interpreted as a convenient representation of the results 
on a limited range of shear rates. In our system, the 
microscopic stress axy is defined by the usual formula 
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(5) 
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The shear rate dependent viscosity is defined by the ratio 
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FIG. 1. Top: Flow curves for various temperatures. 

Full lines are Newtonian behavior (a = 0) and shear-tliinning 
behavior with a = 2/3. Bottom: Flow curves presented in 
the alternative form 77(7), for the same temperatures. 



Figure |^ presents the flow curves 17(7, T) for various 
temperatures. The same data are also represented in 
the alternative form of viscosity versus shear rate. Two 
different regimes can be distinguished in these flgures. 



For T > Tc, the behavior becomes Newtonian at low 
shear rates. 



where 



a ~ 7?o(r)7, 



?7o(r) = lim n(7,r > Tc 

7-^0 



(7) 



(8) 



is the viscosity of the fluid in the linear regime. We 
have checked that it is proppctional to the equilibrium 
relaxation time of the fluidtil, and it is found to di- 
verge as the temperature is lowered as a power law 
VoiT) ^ {T — Tclz^ '*^, in good agreement with equilib- 
rium simulationso. Figure ^ shows that in this temper- 
ature regime, all the data can be rescaled onto a master 
curve of the form 



Vo{T) 



[1 + 7/7o(T)] 



2/3 ■ 



(9) 



7o(r) is not a free parameter, but has the form 70 (T) ~ 
(T — y^) -2.45x3/2^ which ensures that the divergence of 
the viscosity at Tc is a power law r]{j,Tc) ~ 7"^/'^. 




FIG. 2. Scaling behavior of the flow curves 77(7) for various 
temperatures above Tc. The full line is a fit to equation (0). 



Equation describes a crossover between a Newto- 
nian behavior at small shear rates and a power law be- 
havior with T] ~ 7'~2/3 higher shear rates. Remark- 
ably, the same behavior, with the same expoiaept 2/3 
were predicted by the mean-fleld theory of Ref.E3. This 
coincidence, however, is not necessarily significant. A 
similar shear-thinning behavior, rj ^ 7 with expo- 
nents in the range a = 0.5 — 1.0 is obtained in many 
different soft systems, suGk-.as suspensions and concen- 
trated polymer solutiongilu3_ Confined polymer layers 
under constant normal loadEH were also found to exhibit 
very similar shear-thinning behavior. In our case, the 
data could probably be fitted with a slightly different 
value for the exponent, arid, the choice of 2/3 is moti- 
vated by theoretical resultsllj. 
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FIG. 3. Evolution of the effective shear-thinning exponent 
a{T) and of the stretching exponent I3{T). The horizontal ar- 
row shows the value a — 2/3 obtained above Tc. The vertical 
dashed hne is at Tc = 0.435. 



Although our results are obtained on 3 decades of shear 
rates, we cannot report a definite functional form for the 
flow curves. An alternative view of the low temperature 
flow curves is that the effective shear-thinning exponent 
will saturate to the value a(T) = 1 in the (numerically 
unreachable) very low shear rate limit, indicating the ex- 
istence of a yield stress cro(T) = lim^^o o'(7, T) for the 
system at finite temperatures. Hence, we also tried to 
use Eq. (|^), in order to extract a possible value of cro{T). 
Note that this is made possible by the fact that Eq. (Q) 
has more free parameters than a single power law. Fits 
are thus satisfactory, and the values of the yield stress 
obtained at different temperatures in this manner are re- 
ported in Fig. ^. This 'jamming phase diagram' is very 
similar to the one proposed Liu and NagelE3, and stud- 
ied experimentally in Ref.ES. It is also an interesting 
way of thinking of the glass traiisition, and as such was 
theoretically investigated in Ref£3. The 'jammed phase' 
would correspond here to the low-cr, low-T part of the 
diagram. 
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FIG. 4. The (cr, T) plane of the jamming phase diagram. 
The dashed curves are the viscosity contour plots, rj — 20, 
30, 50, 100, 200, 500, 1000 and 2000 (from top to bottom). 
The full line represents the value of the yield stress ao{T), 
extrapolated from Eq. (^. Arrow marks the mode-coupling 
temperature Tc. 

At low temperatures, the curves shown in Fig. |l| 
present a crosssover between two shear-thinning regimes. 
For relatively high shear rates, they can be described by 
the same shear-thinning exponent as above Tc, a = 2/3, 
but at low shear rate, the shear-thinning behavior is 
more marked, a > 1. One can define an effective shear- 
thinning exponent from the slope of the curves in Fig. |l| 
{a = dln?7/dln7)). It crosses over from a ~ 2/3 at rel- 
atively high shear rate, to a larger value for the lowest 
shear rate investigated here. The temperature depen- 
dence of the latter is reported in Fig. ^, which shows 
that a{T) saturates to 1 when T ^ 0, which seems to 
indicate that the system has a finite yield stress in this 
limit only. 



V. MICROSCOPIC CORRELATION FUNCTIONS 

The shear rate may be viewed as a new control param- 
eter to access the glassy phase. As discussed earlier, one 
may expect scaling properties of the correlation functions 
as a function of shear rate to be similar to those observed 
when temperature is-decreased. This basic remark was 
already made in Ref.EJ, but we now have a clear theoret- 
ical context in which these results may be understoodlij. 
In this Section, we briefly consider the influence of the 
shear rate on static correlations, and investigate next the 
scaling properties of the time-dependent correlation func- 
tions. 



A. Static correlation functions 

Static correlations may be characterized through the 
structure factors of the fluid, defined as 

(10) 

We show first in Fig. || the structure factor for 7 = 
10~'^, and T = 0.3 in the three spatial directions. It can 
be seen that the flow, somehow surprisingly, intmduces 
no obvious anisotropy in this structural quantityEj. This 
can be explained by the fact that we are working at shear 
rates which are small, 7T0 <§; 1. A structural anisotropy 
would probably be present for larger shear rates, 7T0 <; 1. 

We present next in Fig. || the evolution of S{kz) for dif- 
ferent shear rates at constant temperature T < T^. The 
'glassy' state is thus approached when the shear rate is 
lowered. We shall see below that for these values of the 
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shear rate, the relaxation time increases by nearly two 
orders of magnitude. Only weak changes in the static 
structure factor are observed in the same shear rate win- 
dow. The situation is thus strongly reminiscent of the 
observation made at equilibrium that no. structural sign 
of the glass transition can be observedcBtJ. 



is the density fluctuation at point x and time t, p being 
the average density of the fluid. For simplicity, we will 
limit ourselves to k vectors orthogonal to the flow direc- 
tions, for which the overall motion of the fluid does not 
affect the correlation function directly. 




FIG. 5. Structure factor at fixed temperature, T = 0.3, 



and shear rate, 7 = 
respect to the fiow. 



10" 



for the different orientations with 




FIG. 6. Structure factor at fixed temperature T — 0.3 and 
different shear rates. The relaxation time increases by ~ 2 
orders of magnitude on the same shear rate interval. 
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FIG. 7. Top: Correlation functions for T = 0.5 > Tc and 
different values of the shear rate, 7 = 0, 10^*, 5 ■ 10^*, 10^'^, 
10~^, 5 • 10^^ (from right to left). Bottom: The slow decay of 
the correlation function can be collapsed if the time is rescaled 
by froi(7). The dashed line is a fit to a stretched exponential 
form, with an exponent P — 0.77. 



To discuss the dynamic behavior of the system, we 
concentrate on the self-part of the intermediate scattering 
function, defined by 

Na 



B. a-relaxation and 'time-shear superposition' Ck[t) = ^exp ^ik ■ [rj{t + Iq) — rj{to 



3 = i 



To investigate the dynamics of the fluid under 
shear, we focus on the number density fluctuations 
{Sp{k,t)6p{~k,t')) where 



(12) 



/ 1 ^ 

M^,t) = /-^<5(a:-r,(t))- 



(11) 



We have also computed this function for the minority 
species, B, with very similar results, that will not be dis- 
cussed here. These correlation functions are displayed 
for two different temperatures and several shear rates in 
Figs. ^ and ||. The wavector is fc = 7.47e2, which cor- 
responds to the peak of the structure factor, see Fig. 0. 
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The overall shape is very similar to what is usually found 
in supercooled systems, namely a first 'microscopic' re- 
laxation which is independent of the control parameters, 
followed by a slow approach to a 'plateau' which we 
will describe as '/3-relaxation' in analogy with the mode- 
coupling terminology. The decay beyond the plateau, or 
terminal relaxation, will be described as 'a-relaxation'. 
Obviously, the a-relaxation time is a strongly decreasing 
function of shear ratec2l. 




be seen that such a scaling is indeed nicely obeyed for 
the a-relaxation of the correlation functions. The scaling 
function is correctly described by a stretched exponential, 

i^fc(x)^exp(-a:^(^)), (14) 

with the stretching exponent < (3{T) < 1. The temper- 
ature evolution of P{T) is shown in Fig. ^, which demon- 
strates that the decay becomes less stretched when the 
temperature is lowered. The stretching exponent seems 
to saturate to the value 1 as T decreases to 0. In this 
figure, a close relationship between this stretching of the 
correlation function and the shear-thinning exponent is 
also clearly apparent. This is physically plausible: (3—1 
corresponds to a pure exponential form, i.e. a single 
relaxation time; /3 < 1, on the other hand, is associ- 
ated with a broader distribution of relaxation times. If 
(3 = 1 — a, there is a single relevant time scale, and the 
scaling rj ^ 7"^ ~ tioi shows that this time scale is sim- 
ply given by the inverse of the shear rate, as would be 
obtained by a naive dimensional analysis. This suggests, 
as an interesting corollary, that a shear-thinning expo- 
nent a < 1 implies the existence of a broad distribution 
of time scales. 
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FIG. 8. Top: Correlation functions for T = 0.3 < Tc and 
different values of the shear rate, 7 = 10~*, 3 ■ 10"'*, 10~^, 
3 • 10"^, 10"^, 3 ■ 10"^, 10"^ (from right to left). Bottom: 
The slow decay of the correlation function can be collapsed if 
the time is rescaled by trci(7). The dashed line is a fit to a 
stretched exponential form, with an exponent 0.95. 

It is a theoretical prediction that the slow decay of the 
correlation function may be collapsed by a simple resGal- 
ing of the time, thp_itime-shear superposition property'113. 
As in equilibriumEj, the relaxation time is defined from 
the relation Cfc(trci) = e~^. The scaling form 

'^'W = (iS^) ■ '''' 

where the function Fk{x) is a master function that may 
depend on temperature, is tested in Figs. and |[ It can 



C. Spatial dependence and /3-relaxation 

The wavevector dependence of the number-number 
correlation function is displayed in Fig. |^. As usualcS, 
the overall shape is qualitatively the same for all wavevec- 
tors, the plateau value and the stretching exponent being 
simple functions of the wavevector. 




FIG. 9. Correlation function for T = 0.3, 7 = 10~^, and 
different values of the wavevector regularly spaced between 
k = 0.93 and k = 24.76 (from top to bottom). 

In mode-coupling thcoriesB'El, an interesting and eas- 
ily testable property of the correlation function is the 
so-called 'factorization property' which characterizes the 
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spatial dependence in the /3-relaxation region. This fac- 
torisation property can be written, for the correlation 
function under study, in the form Ck{t) ^ Qk + hkf{t), 
where qk is the wavcvector-dependent Edwards- Anderson 
parameter, and f{t) a universal function, i.e. indepen- 
dent of the chosen correlation function. More generally, 
it implies that in the /3-regime the ratio 



V{t) - v(t") 



(15) 



where (p{t) is any slow observable [e.g. Cfc(t)], and t, i' 
and t" are in thfij^S-relaxation window, is independent of 
the observableH^EEl. 

This prediction is tested in Fig. which shows that 
the ratios (|l5|) obtained for different wavevectors are in- 
deed independent of the wavevector. Again, thesei data 
are very similar to wiiat is observed at equilibriumEa, and 
in the aging regimecj. 




FIG. 10. The correlation functions of Fig. |^ are rescaled 
in the /3-regime using the mode-coupling-called 'factorization 
property', Eq. @. Times are t' = 0.57 and t" = 14, as 
indicated by the arrows. 



VI. EFFECTIVE TEMPERATURES 

A. Existence of an effective temperature 

Nonequilibrium mean-field or mode-coupling theories 
have shown that whereas correlation functions were the 
only relevant quantities to be studied at equilibrium, it 
is of primary interest to separaihj-study susceptibilities 
in out of equilibrium situationsEj'E3L3. At equilibrium, 
these two families of dynamic functions are related by 
fluctuation-dissipation theorems (FDT). In nonequilib- 
rium situations, FDT is a priori not satisfied. However, it 
now is clearpthat the study of the corrections to the FDT 
are relevantQ. Recent theories of nonequilibrium systems 
have indeed introduced the notion of effective tempera- 
ture associated wit h tJaet modification of the fluctuation- 
dissipation theoremtaEa. This effective temperature is 



defined by the ratio of the response to correlation func- 
tion in the following way. Consider two physical observ- 
ables 0{t) and 0'{t). The correlation function Coo'{t) 
is defined as 

Coo'it) = {0{t + to)0'{to)) - (O(to))(0'(to)), (16) 

where to is an arbitrary initial time and (• • • ) indicates 
an average over different initial times tg. The conjugated 
response function is defined by 



Roo'it) 



S{0{t + to)) 
Sho'ito) 



(17) 



where ho' is the field thermodynamically conjugated to 
0'{t). In a system at thermodynamic equilibrium, the 
FDT reads 



Roo' (t) 



1 dC'oo'jt) 
T At 



(18) 



The more easily accessible physical quantity is the sus- 
ceptibility 



XOO' (t) 



dt'Roo'{t') 



(19) 



which is obtained by applying a small, constant field ho' 
in the time interval [0, i]. In the linear response regime, 
ho' 0, one gets 



Xoo'{t) 



m)-om 

ho' 



(20) 



The equilibrium FDT thus implies a linear relation be- 
tween the susceptibility and the correlation, namely 

XOO' (t) = ^ {Coo- (0) - Coo- it)) . (21) 



In a system out of equilibrium, the FDT is not ex- 
pected to hold. However, in a system invariant under 
time translation, a modified version of Eq. ( |l8| ) can be 

7oo') thro 

dCoo'it) 



j) car 

used to define the functions T^'-' {Coc) througto 



Roo'{t) 



1 



KTiCoo'Y 



dt 



(22) 



The (x) are a priori arbitrary functions of their ar- 
gument, and may depend on which observables 0{t) and 
0'{t) are under study. However, they can be measured 
by following the same linear response procedure as in the 
equilibrium case, so that 



dCoo'it') 



off 
Coo'(O) 



dt' 



(23) 
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The existence of an effective temperature is thus demon- 
strated if a straight line is obtained in a susceptibility- 
correlation plot parameterized by the time. As was 
shown in several casesE^, the same system may exhibit 
different temperatures on different time scales. If this is 
the case, the parametric plot will consist in several dif- 
ferent straight lines, the slope of each line deteijrn.in.ing 
an effective temperature. Cases are also knownOCJ in 
which the parametric plot is ill-defined at short times 
(e.g. because of an oscillatory behavior of the correlation 
function) but yields a well-defined effective temperature 
on longer time scales. It is therefore generally incorrect 
to define an effective temperature through the approxi- 
matcuexpression Xoo'it ^ oo)/Coo'(0), as is sometimes 
doneE3. This corresponds to taking the average slope of 
the parametric plot as an effective temperature, there- 
fore losing all information on the time scale dependence 
of the effective temperature. 

Obviously, the introduction of an effective tempera- 
ture is of 'thermodynamic' interest only if this quantity 
is actually independent of the observables 0{t) and 0{t') 
under considera^oa—Tp'^'^ = Toff- This is true at the 
mean-field leveOEjCJ, and this question will be investi- 
gated in detail below. 

As a first investigation of the existence and behavior of 
the effective temperature and fluctuation-dissipation re- 
lation, we consider the case of single particle density fluc- 
tuations, for which the corresponding correlation func- 
tion has been studied in detail above. Corresponding 
observables are in this case 



^ Na 

0(t) = — ^ Sj exp {ik ■ rj (t)) , 



and 



Na 



O'(0 = 2^£j cos (k-rjit)), 



(24) 



(25) 



where = ±1 is a bimodal random variable of mean 0. 
The relation 



Cfc(i) = Coo'(0, 



(26) 



where the horizontal line means an average over the re- 
alizations of {sj}, easily follows. To compute the suscep- 
tibility, a force 



Fj{k,t) 



d 



ho'O'it) 



(27) 



is exerted on each v4-particle. The response function 
Xkit) [see Eq. (pO|)] is computed by averaging the re- 
sponse obtained in different realizations of the {sj}. 

The time-dependence of the dynamic functions Xk{t) 
and [1 — Ck{t)]/T is displayed in Fig. In that case, 
200 realizations of the {sj} have been considered. This 
figure shows that the equilibrium FDT is obeyed at short 



times, and violated at larger times. From these functions, 
a parametric plot is built, as shown in the inset of Fig. 
It can be seen that the parametric plot can be described 
to an excellent approximation by two straight lines, of 
slopes —1/T (large value of the correlation corresponding 
to short times) and —1/Tcs (small value of the correlation 
corresponding to large times). 

This observation is thus in perfect agreement with the 
behavior predicted theoreticallytj, and was the main re- 
sult in Ref.B. 




FIG. 11. Susceptibility Xfe(i) (squares) and [1 - Ck]{t)]/T 
(circles) versus time, for k — 7A7e^, T = 0.3, and 7 = 10"'^. 
The two curves are superposed when the equilibrium FDT is 
satisfied. Inset: Parametric susceptibility versus correlation 
plot for the same data. The dashed line is a linear fit to the 
small-C part of the data, with Tf.s = 0.65. The full line is the 
equilibrium FDT. 



B. Effective temperature in the plane (7, T) 

The effective temperature is defined as (minus) the 
inverse slope of the parametric susceptibility-correlation 
plot for long times. This effective temperature will, gen- 
erally speaking, depend on how strongly the system is 
externally driven. One may 'reasonably' expect that 
a stronger drive implies a higher effective temperature. 
[See the discussion of the word 'reasonably' at the end of 
the paper.] That this is indeed the case is demonstrated 
in Fig. |l^, where parametric susceptibility-correlation 
plots are displayed for several different shear rates. 

Two different situations are considered in this figure. 
At high temperatures, T > Tc, the effective temperature 
reduces to the bath temperature whea the shear rate is 
decreased. As predicted theoreticallyllj, we find here that 
the effective temperature of the slow modes is equal to 
the bath temperature as long as the system is in its New- 
tonian regime, rj const. This is physically reasonable, 
since the linear regime of rheology corresponds roughly 
to shear rates suctujhat jt^ci < 1, for which the dynamics 
is weakly affectedl30. 
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In the 'glassy' low temperature region, on the other 
hand, the effective temperature also decreases with de- 
creasing shear rate, but saturates for the lower shear rates 
that can be investigated in the simulation at a value dif- 
ferent -Erom the bath temperature, as theoretically pre- 
dictedt3. The existence of such a limiting value is phys- 
ically interpreted by the fact that at low temperature, a 
Newtonian (77 ~ const) or near-equilibrium {Tcs T) 
regime defined as above by the condition 7^101 < 1 does 
not exist, since iici ~ 00. 

The limiting value Teg (7 0) corresponds to the effec- 
tive temperature that was observed in the correspondiag 
aging studies of the same system at this temperatureE3. 
In Ref.E3, the limiting value T^s ^ 0.62 is reported (how- 
ever with much less accuracy than in the present pa- 
per), while we find T^s ~ 0.65 at the same tempexature 
T = 0.3. This is again expected from the theoryll3, and 
confirms the interpretation of Tcff as a signature of the 
geometry of phase space explored by the system on the 
simulation time scale, _i±id£pfindently of the exact proce- 
dure which is foUowectj'EjO. 
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FIG. 12. Parametric plots for T = 0.3 (top) and T = 0.5 
(bottom), and various shear rates. In both figures, the full 
line is the FDT, and has a slope — l/T. The dashed lines are 
linear fits to the small-C part of the data, for 7 > 0. The 
wavevector is fc = TAYe^ for all the curves of this figure. 



VII. OBSERVABLE INDEPENDENCE OF T^a 

As we already mentioned, the notion of an effective 
temperature derived from the fluctuation-dissipation re- 
lation is much more relevant if it can be shown not to 
depend on the observable under consideration. That Tes 
is observable independent is indeed on ci o f the crucial 
predictions of mean-field approacheaialiSiHa. In this Sec- 
tion, we explore this important issue by computing the 
susceptibility-correlation parametric plots for several dif- 
ferent observables. Since the computation of response 
functions is numerically very demanding, the method we 
used consisted in computing very accurately the effective 
temperature from the scattering function for k — lAlez- 
^From the data in Fig. [ill, one has the estimation Tcff ~ 
0.65 for T = 0.3 and 7 — 10~^. We can then compute 
for the same parameters (T, 7) the response associated 
with other observables, and check whether these data 
are compatible with this value of Teff- Since we have ex- 
plained in some detail the generic numerical procedure 
to measure T^g in the previous Section, less details will 
be given here. Instead, we carefully discuss the possi- 
ble experimental realizations of our measurements. Note 
that in experiments, a simultaneous measurement of the 
response and of the correlation is v£f¥-difficult, and has 
been achieved only in very few caseafjcZl. While usual di- 
electric, magnetic or mechanical methods yield response 
functionS|r|S|Cattering experiments probe only correlation 
functionsEra. 



A. Position correlations for different wavevectors 

As a first check, we have investigated the wavevector 
dependence of the parametric plot for the self-part of the 
intermediate scattering function. The results, displayed 
in Fig. |l^, are indeed perfectly consistent with a wavevec- 
tor independent effective temperature. Parametric plots 
for different wavevectors differ only through the value 
of the correlation function at which the crossover from 
bath to effective temperature takes place. This value 
marks the crossover between 'short', equilibrated and 
'long', out of equilibrium time scales. It also corresponds 
to the plateau value qk in the correlation functions dis- 
played in Fig. 0. Note that the wavevectors studied here 
cover a range of length scales roughly between O.baAA 
and 'iaAA- That these different length scales have the 
same off-equilibrium behavior, as observed through the 
effective temperature, is not a priori obvious. It is in 
opposition, for instance, to the simple view of the sys- 
tem rapidly equilibrating its small length scales, while its 
large length scales need much more time to reach equilib- 
rium, a generic schpne which is put forward in systems 
such as spin glassesE3. 

The 'chemical' dependence of the effective temperature 
can also be investigated by computing the correlation 
function Ck{t) and the corresponding response, associ- 
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ated with A and B particles of the Lennard- Jones binary 
mixture. The results, displayed in Fig. again confirm 
that the slope of the parametric plots does not depend 
on the particular chosen observable, all the plots being 
well compatible with the same Teff = 0.65. 




0.2 0.4 0.6 0.8 



FIG. 13. Comparison of the parametric plot for T — 0.3 
and 7 = 10~^, and different wavevectors in the direction 
z. The full line is the equilibrium FDT of slope -1/0.3, the 
dashed lines have slope -1/0.65. 
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FIG. 14. Comparison of the parametric plot for T — 0.3 
and 7 = 10"'^, for the two types {A and B) of Lennard- Jones 
particles. The data are compatible with the same Tcff al- 
though the statistics for B particles, the minority species, is 
slightly not as good as for A particles. The full line is the 
equilibrium FDT of slope -1/0.3, the dashed lines have slope 
-1/0.65. 

Finally, the collective correlation and response func- 
tions can also be computed — albeit with a bit less ac- 
curacy than the incoherent ones. This correspond to ob- 
servables (^) and ( p5| ) but taking Ej = 1 for all j. The 
resulting parametric plots are shown in Fig. |l^, for two 
diff'erent wavevectors. Again, slopes are perfectly com- 
patible with equilibrium FDT for short times, and with 



Toff = 0.65 for longer times. 

It is interesting to consider how such response- 
correlation plots could be obtained experimentally in soft 
condensed matter systems. While the correlation func- 
tions are reasonabbfraaaily obtainable through light scat- 
tering experimentsOijO, the same is not true of response 
functions. In order to obtain the latter, one has to be 
able to manipulate the particles through some externally 
applied potential, modulated at the same wavevector as 
used in the light scattering experiment. One suggestion 
here would be to use some non-index matched tracer 
particles in an index matched colloidal suspension. The 
tracer particles would then be sensitive to the intensity 
of the local electric field. The same efi'ect is used, e.g., 
in optical tweezers. An interference pattern would actu- 
ally realize the modulated external potential considered 
in this Section. Reading of the response could then be 
obtained from the scattering at the wavelength corre- 
sponding to this pattern. 
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FIG. 15. Parametric plots from the coherent density fluctu- 
ations for two wavevectors k — 7A7e^ and k = 12.62e2 . The 
full line is the equilibrium FDT of slope -1/0.3, the dashed 
lines have slope -1/0.65. 



B. Einstein relation for the Lennard- Jones particles 

A quantity which is not directly related to the incoher- 
ent correlation function is the mean square displacement 
of a tagged particle. The response function that is asso- 
ciated with this quantity is the displacement induced by 
applying a small, constant external force to this tagged 
particle. Both quantities are linked by a fluctuation- 
dissipation theorem, called the Einstein relation in that 
context. To increase the statistics, we compute the mean 
square displacement of all particles of the same type, for 
example for the particle A, we measure 

^ ^ E ([^^-(^ + ^o) - ^^ (^o)] 7 ■ (28) 
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Again, we focus on the direction transverse to the flow, 
for computational simplicity. The conjugated response 
function is then 



Xo{t) 



1 



NaFo 



Na 



to) 



(29) 



where Zj (t) is the transverse position of particle j at time 
t. The force Fj = SjFq is applied in the direction z 
between times to and I+Iq. The FDT reads then xo(0 = 

miT. 




FIG. 16. Mean square displacement ( BSf) n ormalized by the 
temperature, and induced displacement (p9|) as functions of 
time. At equilibrium, the FDT implies the equality of both 
quantities. Full lines are the power law behavior of ballistic 
{t^) and diffusive (t) regimes. 



At short times, the mean square displacement is ballis- 
tic, A{t) ^ . This ballistic regime is, as usual, followed 
by a diffusive regime A(t) t. K more extensive dis- 
cussion-of this quantity in the shear flow may be found 
in Ref.Ell. The time behavior is similar for the induced 
displacement, as shown in Fig. |l^. Hence, by construc- 
tion, a linear relation is obtained in the diffusive regime 
in the parametric plot of xo versus A. We confirm in 
Fig. ^ that, again, the effective temperature defined in 
this way is compatible with the value Teg = 0.65 already 
obtained above. This is true for the two types of the 
Lennard-Jones mixture, see Fig. |l^. 

One may therefore define T^fi by simply using the re- 
lation 

Teff = DC (30) 

between the diffusion constant D and the friction coeffi- 
cient C which are defined, as usually, as 

1 



(31) 



for times which are in the_diffusive regime. This method 
was already used in RefsEjo to extract-an effective tem- 
perature in aging systems, and in Ref.H in the context 
of sheared granular materials. 



Its interest, especially in view of experimental realiza- 
tions of these measurements, lies in the fact that the full 
time dependence of the correlation and response func- 
tions is not needed to extract the effective temperature. 
Again, experiments could be considered if tracer particle 
sensitive to an external force field (e.g. magnetic parti- 
cles) could be introduced into the system, and their mo- 
bility ( measured together with their diffusion constant 
D. 
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FIG. 17. Parametric plot for T = 0.3 and 7 = 10"^, for 
the self-diffusion of the two types of Lennard-Jones particles. 
We have represented the quantity xo versus a; — A, where 
X — 0.6 for B particles, x = 0.8 for A particles, for graphical 
convenience. These curves are compared to the data obtained 
for the self part of the intermediate scattering function at 
fc — 7.47ej. The full lines are the equilibrium FDT of slope 
-1/0.3, the dashed lines have slope -1/0.65. 



C. Stress fluctuations. 

A completely different observable, of relevance to all 
flow situations, is the stress defined in Eq. (||). In this 
Section, we study thus the case in which the two observ- 
ables 0{t) and 0'{t) are equal to the diagonal stress in 
the direction perpendicular to the velocity and velocity 
gradient, 0{t) = 0'{t) = (Tzz{t)- To add a field con- 
jugated to the normal stress, a compression 5Lz of the 
simulation box is realized by rescaling all particle coor- 
dinates at time t = Q. In the limit of small compression, 
5Lz/Lz <C 1, the Hamiltonian H of the system becomes 
indeed 



H{5Lz) = HiSLz 



0) + —^Vdzz 




(32) 



where V = L^LyL 



is the volume of the system. This 
proves that dL^L^Ly is the field thermodynamically con- 
jugated to ffzz and allows us to access another check of 
the observable independence of Teff. 
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The resulting parametric plot is shown in Fig. Al- 
though the data are more noisy than in the density cor- 
relation case, the FDT violations are very similar, and 
the value of Toff is well compatible with all the others. 

Again, it may be possible to carry out the experiment 
corresponding to this simulation. Experimentally, the 
off-diagonal component of the stress (rather than a diag- 
onal one) would be used as the observable. In fact, pre- 
liminary results in this direction, using an extremely sen- 
sitive rheometer, have been obtained on aging colloidal 
systems by Bellon and Cilibertdlj. The superposition of 
a continuum flow in order to obtain the fluctuations in 
a sheared system, however, would represent a difficult 
experimental challenge. 
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FIG. 18. Parametric plot for T = 0.3 and 7 = 10"^, for 
the stress fluctuations. The full line is the equilibrium FDT 
of slope -1/0.3, the dashed line has slope -1/0.65. 



VIII. FURTHER EXPERIMENTAL PROPOSALS: 
EQUIPARTITION THEOREM FOR THE SLOW 
MODES 

As mentioned above, tracer particles can be used to 
probe the dynamics of the fluid. In particular, the Ein- 
stein relation between mobility and diffusion leads easily 
to the definition of an effective temperature, which ap- 
pears to be identical to the one obtained from various 
susceptibility-correlation functions. 

A different use of tracer particles is considered in this 
Section, in which we investigate the dynamics of tracer 
particles with a mass large compared to that of the fluid 
particles. The idea is the following. The definition of Toff 
as a ratio between correlation and response function has 
been shown to imply that, if the fluid is used as a ther- 
mal bath to equilibrate a subsystem (or 'thermometer') 
of typical time scale tg ~ t^-d, then the Boltzmann weight 
for the subsystem will be givep-by exp(— _Bs/Tcff), with 
Eg the energy of the subsystemE^. The thermometer does 
not measures the microscopic temperature, but the effec- 
tive temperature associated with its characteristic time 



scale. In that sense, Tcs can correctly be described as an 
'effective temperature'. 

Here, we propose a concrete realization of this situation 
in which the thermometer is in fact a heavy particle, with 
a mass mtr 3> m. Tuning the mass of the particle allows 
to control its Einstein frequency uje, defined byo 



0JE^^^Jd3rg{r)VW{r), (33) 

which essentially characterizes the frequency of vibration 
for the particle in the cage formed by its neighbors. For 
particles of type A with mtr = 1, ujetq ~ 20 at T = 0.5. 

In our simulations, we have considered 10 massive 
tracer particles with masses in the range mtr G [Ij 10^]. 
The tracer particles are otherwise identical to A parti- 
cles, so that the equilibrium structure of the fluid is the 
same as with normal, 'light' particles. Since the oscil- 
lation frequency is inversely proportional to the square 
root of the mass, the heaviest particles will have an Ein- 
stein frequency typically thousand times smaller than the 
light ones, which puts their oscillation period in the range 
of typical relaxation times for the sheared supercooled 
fluid. Note that the Einstein frequency is entirely deter- 
mined by static properties, and always determines the 
short time behavior of the velocity autocorrelation of the 
particle. Therefore it is a relevant quantity also for the 
description of heavy particles, moving in a fast environ- 
ment. 
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FIG. 19. Time dependence of the quantity mtr«z in par- 
ticular trajectories of tracer particle of mass mtr ~ 10^, 
10^ • • • 10® (from bottom to top). Curves have been vertically 
shifted for clarity. 

Fig. |l^ compares the typical trajectories for heavy and 
light particles. The heavy particles are indeed seen to 
have a much lower oscillation frequency, and therefore 
their dynamics can be seen as involving a low frequency 
filter coupled to the fluctuations of the host fluid, which 
corresponds t0j_the slow, or 'effective thermometer' con- 
sidered in Rcf.tj. In order to get the temperature from 
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this thermometer, we compute the average mean square 
velocity in the direction transverse to the flow for these 
tracer particles. The results are shown in Fig. From 
this figure, it is clear that light particles have a mean 
squared velocity given by the usual bath temperature 
(T — 0.3 in this case), which crosses over to the effec- 
tive temperature Tcs — 0.65 for the heaviest particles, 
mtr ~ 10^ — 10^. Intermediate masses, rritr ~ 10'* — 10^ 
corresponds to times scales located in the plateau re- 
gion of the correlation which marks the crossover between 
equilibrium and off-equilibrium behaviors. This crossover 
is usually entirely hidden in the 'breaking points' near the 
value Cfc ~ Qk in previous parametric plots. 
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FIG. 20. Mass dependence of the mean kinetic energy in 
the z direction for T = 0.3 and 7 = 10~^. Horizontal lines 
are T — 0.3 and Tag = 0.65. Error bars are evaluated from 
tracer to tracer fluctuations. 




FIG. 21. Probability distribution function of the velocity 
Vz for various massive tracers, as indicating in the Figure. The 
full lines are fits to the form (^) demonstrating the Gaussian 
shape for all the values of the mass investigated. Values of TcS 
are Tcff = 0.3 (equilibrium, mtr = 10'^), 0.44, 0.54 (crossover, 
mtr = 10* — 10^) and 0.65 (asymptotic effective temperature, 
mtr = lO'' - 10'^). 



This implies that a generalized equipartition theorem 
holds for the heaviest tracers, with the effective temper- 
ature replacing the usual temperature, i.e. 



1 



1 



(34) 



Moreover, we show in Fig. ^ that P{vz), the probability 
distribution of the velocity, is well approximated by a 
Gaussian shape. 



■ exp 



(35) 



This in turn implies that the usual Maxwellian shape 
of the velocity distribution is recovered, again with Tcff 
replacing T. 

The result, although expected on general grounds, is 
physically quite surprising. Again, it could be tested 
against experiments involving for instance colloidal parti- 
cles, the tracers, in a complex fluid, e.g. polymeric. The 
important factor here is not the size of the tracer parti- 
cle, but rather its mass that should be much larger than 
that of the fluid constituantsE^. One could also investi- 
gate the rotational degrees of freedom of non-spherical, 
'slow' tracer particles. 

This result opens the way for new, simple determi- 
nations of effective temperatures in out of equilibrium 
glassy materials. Indeed, it has to be noted here that no 
'complex' dynamic functions such as correlations or sus- 
ceptibilities are needed. This suggests that it would be. 
extremely interesting to reproduce Perrin's expcrimcntL£l 
on barometric equilibrium in a sheared colloidal suspen- 
sion. Although we have not tested numerically this situa- 
tion, we expect indeed that the barometric equilibrium of 
heavy particles inside a horizontally sheared fluid should 
also be ruled by the effective temperature of the fluid, 
higher than the room temperature. 



IX. DISCUSSION 

In this paper, we have investigated the correlation and 
response functions for a supercooled fluid undergoing 
steady shear flow. The ingredients essential to our study 
are the separation of time scales existing in the fluid at 
rest between the slow a-relaxation and the microscopic 
times, and the energy input at large scales provided by 
the shear flow. With these minimal ingredients, it is pos- 
sible to carry out a meaningful comparison with theo- 
retical results obtained for this noncquilibrium dynamics 
at the mean-fleld, or schematic mode-coupling level. We 
also believe that the results should be of general rele- 
vance to many complex systems undergoing shear flow, 
for which the same ingredients are present. 

We have first studied the nonlinear rheology of the 
system, which is indeed quite similar to what is observed 
in many complex systems. The system exhibits shear- 
thinning for shear rates exceeding the a-relaxation time. 
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The results can be rescaled on a common curve as long 
as temperature remains higher than the critical mode- 
coupling temperature for this system. Below Tc, the 
shear-thinning behavior may alternatively be interpreted 
in terms of a temperature-dependent shear-thinning ex- 
ponent, or a temperature-dependent yield stress. The 
fact that tr is nearly independant of the shear rate at low 
temperatures and shear rates makes the Lennard-Jones 
mixture a good candidate to study possible mechanical 
instabilities in Bingham fluids at a microscopic level. 

In the shear flow, the system dynamics becomes sta- 
tionary. This allows an accurate determination of the 
time dependent correlation functions, as opposed to the 
case of aging systems, for which this determination is 
more difficult. We showed that several generic features 
of correlation functions observed in supercooled systems 
in the vicinity of Tc are stilL^resent out of equilibrium, as 
was predicted theoreticallyll3. Among those, the two-step 
(/3 and a) relaxation, the stretching of the a-relaxation 
and the factorization property which characterizes the 
spatial dependence in the /3-relaxation regime. 




0.2 0.4 0.6 0.8 1 

Correlation 



FIG. 22. Summary of our results for the effective tempera- 
ture. We show here the 14 different susceptibility-correlation 
parametric plots described in the previous sections. They are 
all consistent with the same effective temperature for the slow 
modes. The full lines are the equilibrium FDT of slope -1/0.3, 
the dashed line have slope -1/0.65. 

Finally, we have investigated in great detail the vi- 
olations of the fluctuation-dissipation theorem in this 
nonequilibrium system. Our results, which are summa- 
rized in a spectacular way in Fig. |2^, are consistent with 
the expectation -based both on theoretical arguments 
and on earlier results from aging studies- that an ef- 
fective temperature, independent of the observable un- 
der consideration, can be defined for the slow modes of 
the fluid. For several of the observables we considered, 
we suggested experimental probes that could be used to 
measure similar quantities in complex fluids. 

We also presented a numerical realization of an 'effec- 
tive thermometer' that probes this effective temperature, 



using heavy tracer particles. The reading of the value of 
the effective thermometer follows from the generalization 
of the equipartition theorem to this nonequilibrium situa- 
tion. This suggest a simple and elegant method to access 
the effective temperature. Again, we suggested experi- 
mental protocols that could realize these measurements. 

It has to be noted that the notion of an effective (or Ac- 
tive) temperature is a long-standing issue in the context 
of glassy materials (see e.g. Ref.El). However, our results 
clearly show that the definition Eq. (|l|) from a nonequi- 
librium fluctuation-dissipation has all the-propcrtics that 
has to be captured by such a quantityll3. (1) It is the 
temperature measured by a thermometer, and it includes 
naturally the separation of time-scale of glassy systems, 
as shown by our 'tracer experiment'. (2) It can be mea- 
sured in a purely dynamic way, without making reference 
to any equilibrated, or extrapolated, state. (3) It is in- 
dependent of the chosen observable, as we have clearly 
demonstrated. None of these points is satisfied by previ- 
ous definitions of a fictive temperatureEJ. Note that Toff 
also captures the essential phenomenological idea that 
when a system is sheared more vigorously its effective 
temperature increases. We thus emphasize here the huge 
physical interest of this quantity, well beyond the simple 
check of a sophisticated theoretical prediction. 

All these results give strong support to the theoretical 
scenario elaborated from schematic mode-coupling the.- 
ories to describe the rheology of soft glassy materialsE^. 
This shows that although these theories give only a rather 
crude description of a sheared fluid, they also capture its 
essential features. It is important to note that some of 
these features are present in no other theory, such as the 
emphasis put on the study of response functions, and the 
resulting FDT violations. 

At the theoretical level, the next step to be performed 
is of course to give a more 'microscopic' derivation of the 
nonequilibrium mode-coupling theory of RcfJ13 for the 
study of glassy materials under shear. In this respect, 
the results obtained by LatzH in its work on the aging 
regime are a very important step. 

However, we believe that it now is time to perform sys- 
tematic experiments to test quantitatively existing theo- 
ries of the nonequilibrium dynamics of glassy materials. 
In that sense, the situation is more or less similar to the 
mid-eighties when schematic equilibrium mode-coupling 
were already derived, but with little experimental con- 
firmation of its main features. This is why we tried, as 
much as possible, to suggest experimental counterparts 
to our numerical measurements. We hope that our find- 
ings and suggestions in the present paper will motivate 
further experimental work in the field. 
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